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Enhanced  Doppler  Picture  Velocimetry  (DPV)  for  Planar 
Velocity  Measurements  in  High  Speed  Shock  Tunnel  Flow 
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Abstract:  A  technique  for  visualizing  a  velocity  field  in  an  entire  plane  has  been  developed  by 
taking  so-called  “Doppler  Pictures”  using  Michelson  interferometry.  With  the  Doppler  Picture 
Velocimetry  (DPV)  information  about  the  instantaneous  and  local  velocities  of  tracers  passing 
through  a  light  sheet  are  available.  The  tracer  particles  are  illuminated  by  a  laser  light  source  and 
the  frequency  of  the  scattered  light,  in  case  of  moving  particles,  is  shifted  by  the  Doppler  effect. 

This  small  Doppler  shift  of  the  frequency  of  the  light  scattered  by  tracer  particles  is  transformed  by 
a  Michelson  interferometer  into  varying  light  intensities  as  interferometer  output.  Therefore,  the 
light  intensity  distribution  on  the  Doppler  picture  gives  information  on  the  Doppler  frequency  shift 
and  in  consequence  on  the  speed  of  the  tracer  particles  crossing  the  light  sheet  plane.  The 
technique  for  taking  and  processing  the  Doppler  picture  images  was  enhanced  in  the  last  years  and 
the  status  of  progress  of  the  DPV  method  will  be  described  first  with  the  tests  using  a  rotating  disc 
and  a  free  jet  for  DPV  technique  calibration,  and  secondly  with  an  application  in  high  speed  fluid 
mechanics  showing  the  velocity  distribution  in  a  light  sheet  plane  crossing  a  supersonic  wedge  flow 
generated  in  ISL’s  high  energy  shock  tunnel  STB.  The  Doppler  picture,  taken  with  the  enhanced 
set-up  by  a  CCD  camera,  is  pixel-wise  computer-processed  for  velocity  presentation  using  an 
image  processing  software  which  is  specially  developed  at  ISL  for  DPV  purposes. 

Keywords:  Light  sheet,  Tracer,  Velocity  visualization,  Michelson  interferometer,  Doppler  picture 


1.  Introduction 

Some  measuring  techniques  have  been  developed  for  taking  a  picture  of  the  velocity  distribution,  for  example  in 
a  light  sheet  crossing  a  gas  flow.  A  well-known  technique  is  the  Particle  Image  Velocimetry  (PIV),  which  is  an 
often  used  tool  for  this  purpose  (Merzkirch,  1990).  Also  Doppler  Global  Velocimetry  (DGV)  is  known,  using  an 
amplitude  light  filter  for  analysing  the  Doppler  shift  of  the  frequency  of  the  light  scattered  by  tracer  particles 
(Meyers,  1992).  At  ISL  a  special  Michelson  interferometer  was  designed  (Oertel  et  al.,  1981)  to  transform  the 
Doppler  frequency  shift  of  the  light  scattered  by  tracer  particles  passing  a  light  sheet  into  a  shift  of  the  light 
intensity  leaving  the  interferometer.  The  changes  in  the  light  intensity  distribution  can  be  directly  related  to  the 
velocities  of  the  tracer  particles.  This  method,  called  Doppler  Picture  Velocimetry  (DPV),  was  further  developed 
during  the  last  years  (Seiler  et  al.  ,1987,  1991,  1998)  and  the  actual  status  of  the  DPV  set-up  and  the  technique 
for  processing  the  interference  image  of  the  Doppler  picture  will  be  described  herein.  Some  details  on  the 
progress  in  image  processing  were  already  described  by  Seiler  et  al.  (1999)  and  Leopold  et  al.  (2000). 


2.  Principles  of  DPV  technique 

2.1  Doppler  effect 

In  order  to  measure  the  tracer  velocity  v ,  the  well-known  relation  of  single-beam  velocimetry  is  used  as  given  in 
Equation  (1)  with  the  unit  vectors  eL  and  eD  .  In  Figure  1  the  vector  eL  denotes  the  direction  of  illumination  by 
the  light  source  (L),  and  the  vector  eD  is  directed  from  the  tracers  to  the  Michelson  interferometer  detector  (D). 
0  is  the  angle  between  the  two  unit  vectors  eL  and  eD .  The  velocity  component  of  the  flow  vector  v  is 
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measured  in  the  direction  of  the  difference  vector  eD  -eL  which  forms 
the  angle  7  to  the  velocity  vector  v ,  as  shown  in  Figure  1.  The  Doppler 
frequency  shift  can  be  deduced  with  the  Doppler  effect  relations  as 
follows: 


-i  v  •  9 
2  —  cos  7  am 


(1) 


The  laser  light  source  frequency  is  denoted  by  vL.  By  Doppler  effect 
that  frequency  is  shifted  to  Vp  and  the  Doppler  shift  dvL  =  Vd  -  Vl  is 
analyzed  via  the  Michelson  interferometer. 


2.2  Michelson  interferometer 

A  plane  I  crossing  the  gas  flow  is  illuminated  by  a  light  sheet  with  frequency  vL.  In  order  to  visualise  the  velocity 
distribution  in  this  plane  X,  tracer  particles  are  seeded  into  the  flow  and  the  light  scattered  from  these  particles  is 
focused  into  a  specially  designed  Michelson  interferometer  (Seiler  and  Oertel,  1983),  see  Figure  2. 


Fig.  2  Whole  field  Michelson  interferometer 


The  Michelson  interferometer  in  Figure  2  consists  of  a  beam  splitter  cube  T  with  a  50%  reflecting  mirror,  the 
two  mirrors  Si  and  S2  and  the  glass  block  G  between  S2  and  T.  The  scattered  light  coming  from  the  object  plane 
Z  and  passing  through  the  lens  L,  is  divided  by  the  beam  splitter  T  into  two  parts  of  equal  intensity.  These  two 
parts  are  focussed  by  the  lens  on  S!  and  S2,  respectively.  With  the  lens  L2  the  image  of  the  object  plane  X  on 
Sj  and  S2  is  transferred  to  the  image  plane  X'.  The  Michelson  interferometer  works  as  a  spectrometer  for 
detecting  the  Doppler  shift  of  the  frequency  of  the  light  scattered  by  tracer  particles  and  its  mechanism  is 
described  in  the  following.  In  the  image  plane  X'  the  light  coming  from  the  two  interferometer  legs  (1)  and  (2) 
interferes  and  the  light  intensity  I  at  all  image  points  P'  in  X',  as  a  result  of  that  interference,  depends  on  the  phase 
difference  Atp,  see  Figure  2,  originated  in  the  two  legs  of  the  Michelson  interferometer  as  given  in  the  following 
relation,  with  I  the  light  intensity  occurring  by  light  interference  on  the  image  plane  X',  and  I0  the  light  intensity 
as  it  is  scattered  by  the  tracer  particles: 

I=I0  cos2  ~  ,  normalized  relation:  — =cos2  —  with  0<— <1  .  (2) 

2  Io  2  I0 

The  phase  difference  A<p  is  a  function  of  three  quantities:  the  optical  path  difference  A<|>,  with 

A*=2n|l2-l,|,  (3) 

the  light  source  frequency  vL  and  the  speed  of  light  c0  in  vacuum,  as  follows: 
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A(p=v 

2k  l  c0 

In  differentiating  A(p  in  Equation  (4)  to  dvL,  the  expression  (5)  is  obtained.  Equation  (5)  explains  that  frequency 
changes  dvL,  due  to  the  Doppler  effect,  are  directly  proportional  to  variations  d(Acp)  of  the  phase  difference  A<p. 

d(|2)=^dvL.  (5) 

ZK  Cq 

The  light  scattered  from  the  light  sheet  plane  £  is  Doppler  shifted  to  vD  in  the  case  of  tracer  displacement  caused 
by  flow  movement.  The  frequency  changes 

dvL=vD-vL.  (6) 

of  the  scattered  light  coming  from  the  plane  £  are  transformed  into  changes  d(A(p),  see  Equation  (5),  of  the  phase 
difference  Acp  between  the  two  legs  of  the  Michelson  interferometer.  As  already  said,  in  Equations  (1),  (4),  (5) 
and  (6)  the  index  (L)  denotes  the  light  source,  index  (D)  the  detector,  i.e.,  the  Michelson  interferometer.  The 
d(Acp)  change  in  Equation  (5)  causes  variations  dl  of  the  interference  light  intensity  I  in  the  plane  £',  as: 

dI=~^-sinA(p-d(A(p)  .  (7) 

Therefore,  the  light  intensity  distribution  I  in  £'  (Doppler  picture),  relative  to  the  initially  adjusted  I  as  reference 
image,  informs  on  the  frequency  shift  dvL  which  is  proportional  to  the  tracer  speed:  relation  see  Equation  (1).  I.e., 
at  all  points  P’  in  £'  dl  is  given  as:  dl  =  I  (reference  image)  - 1  (Doppler  picture). 

The  Michelson  interferometer  is  adjusted  initially,  before  flow  onset,  to  interference  fringes  by  turning  the 
mirror  Si  by  a  small  angle  (3  as  depicted  in  Figure  2.  It  would  be  preferable  to  adjust  the  interferometer  to  a 
homogeneously  light  intensity  output,  which  is  a  constantly  distributed  light  intensity  I  in  plane  £',  but  is  not 
recommendable  in  practice  because  the  optical  components  such  as  mirrors,  lenses  and  other  parts  used  have  only 
a  limited  optical  accuracy.  Therefore,  the  best  alternative  is  to  adjust  the  Michelson  interferometer  initially  to  the 
easily  obtainable  interference  fringe  pattern  mentioned  (Seiler  and  George,  1 986). 


3.  Tests  for  DPV  system  calibration 

3.1  Rotating  disc 

The  periphery  of  a  rotating  disc  is 

illuminated  by  a  1  W  ArMaser  from  a 

direction  nearly  perpendicular  to  the  axis 

of  rotation.  The  scattered  light  is  observed 

in  the  hack  scattered  direction,  see  motor 

principle  sketch  in  Figure.  3.  cowling 


For  the  rotating  disc,  the  velocity  vector  of 

the  “scattering  surface”  is  tangential  to  the  Ar+-laser 

surface.  This  tangential  velocity  can  be 

divided  in  two  vectors:  one  vector  in  the  measured 

direction  of  observation  (Michelson  velocity 

interferometer)  and  one  perpendicularly  to 

that,  which  causes  no  Doppler  effect. 


visualized  zone 


Michelson  interferometer 


Fig.  3  Set-up  of  the  rotating  disc 


The  measured  velocity  component  v  is,  according  to  Equation  (1),  proportional  to  the  Doppler  shift 
dv  L  =vd  ~vL  which  is  processed  by  the  Michelson  interferometer  for  velocity  determination. 
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In  the  Doppler  picture  technique,  the  interference  fringes  have  been  adjusted  by  turning  the  mirror  Si  by  a 
small  angle  p,  see  optical-set-up  in  Figure  2.  Then  straight  and  parallel  interference  fringes  appear,  see  Figure  4a. 
Two  fringe  light  intensity  patterns  are  visible  which  depend  on  the  phase  difference  Arp  according  to  the  cos2- 
distribution  in  Equation  (2).  On  the  left  hand  side  of  the  Doppler  picture  in  Figure  4a  the  interference  fringe 
system  is  formed  by  the  light  scattered  from  the  motor  cowling  at  rest,  and  on  the  right  the  fringe  system  comes 
from  the  light  scattered  from  the  rotating  disc  periphery.  The  disc  axis  is  at  the  top  of  Figure  4a. 

In  order  to  evaluate  the  velocity  at  each  pixel  of  the  CCD  camera,  at  first  each  of  the  two  interference  fringe 
systems  is  normalized  with  the  actual  light  intensity  distribution  I0  (Figure  4b),  with  I/I0  according  to  Equation 
(2).  I0  is  captured  by  the  camera  CCD1,  see  Figure  5.  Imperfections  in  the  optical  parts,  non-uniformity  of  pixel 

distribution  and  disorientations  make  the  exact  overlay  of  the  interference 
(CCD2)  and  intensity  images  (CCD1)  very  difficult  (Meyers  1992). 

Therefore,  after  each  alignment  of 
the  optical  system,  several  images  of  the 
calibration  card  placed  in  the  light  sheet 
plane  are  recorded.  The  calibration  card 
consists  of  4  different  lines  and  36  grid 
points,  Figure  6.  The  4  lines  were 
detected  with  an  algorithm  developed  at 
disc  ISL  by  Schneider  (1987).  By  means  of 
the  information  concerning  the  position 
and  the  orientation  of  these  lines,  the 
images  can  be  correctly  aligned  to  each 
other. 

Fig.  4a  CCD2  interference  image  Fig.  6  Calibration  card 


Fig.  5  Intensity  camera  CCD1  and  DPV  camera  CCD2 

Several  sets  of  linear  equations  are  established  in  order  to  move  the  grid  points  of  the  calibration  card  from 
their  imaged  to  their  overlayed  positions.  To  get  an  accurate  transformation  of  the  gray  values,  the  interpolation 
technique  described  by  Schumacher  (1996)  is  applied. 
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Because  of  contrast  variations  in  the  DPV  pictures,  Equation  (2)  is  not  really  fulfilled  with  0  <  I/I0  ^  1  over 
the  whole  interference  fringe  picture,  i.e.,  the  maximum  gray  values  as  well  as  the  minimum  values  fluctuate.  The 
largest  G ^  and  smallest  Gmim  gray  values  are  determined  with  every  pixel  on  the  CCD  chip  in  its  (x,  y)-plane 


and  a  second  normalization  technique  is  applied,  using  the  G„ 


G(x,  y)  =  G^  +  (G^  -G^  )cos' 


.  values  in  the  following  way: 


In  using  the  expression 

i(xj)=f4rk=cos2f’  (9) 

^max  'Jmin  ^ 

the  equation  (9)  is  transformed  to  Arp  which  results,  in  case  of  motionless  motor  cowling  as  reference,  to: 

AcpL  =arccos(2l(x,y)L -l)±  n  n=l,2,3,...  (10) 

I(x,  y)L  is  determined  from  Figure  4a,  left  hand  side  disc  fringe  pattern.  The  periodical  fringe  pattern  cause  an 
increase  of  n.  For  the  right  hand  side  of  the  interference  fringe  distribution,  Arp  is  determined  in  the  same  way: 

A(pD  =arccos(2l(x,y)D -l)  ±  n  n  =  l,2,3,...  (11) 

Finally,  assuming  a  small  difference  between  the  laser  frequency  vL  and  the  Doppler  shifted  one  vD,  the 
Doppler  shift  dvL  results  from  Equation  (5)  as: 

dVl  =^dM  ,  (12a) 


d(A<p)  =  A(pD  -  AcpL 

Figure  7  shows  the  Doppler 
picture  in  terms  of  the  phase 
difference  Arp  treated  as  described  by 
Equation  (10)  and  Equation  (11)  as 
pseudo  color  distribution.  Figure  8 
represents  the  velocities  measured  at 
each  pixel  of  the  Doppler  picture 
image  applying  Equation  (1)  and 
Equations  (12a, b).  The  disc  axis  is 
located  at  x  =  23  mm  at  zero  speed 
measured.  There  are  small 
differences  between  theory  and 
Doppler  picture  measurement  which 
increase  towards  the  upper  part  of  the 
disc  for  x  -»  0.  The  deviation  is 
systematic  and  probably  caused  by 
the  inaccuracy  of  the  calculated 
speed  distribution  (=5%)  and  by  the 
error  of  Doppler  picture  processing 
(maximum  error  of  is  10%),  which 
increase  with  the  decreasing  fringe 
contrast  towards  x  0. 


A<Pl 

Eq.(10) 


A(pD  d(A(p) 
Eq.(ll)  Eq.(12b) 


disc  axis 
x  =  23  mm 


- theory 

0  greyvalues 


0  50  100  150  200  250  300 


Fig.  7  Pseudo  color  Doppler  pictures 
(left:  motionless,  right:  rotating  disc) 


Fig.  8  Velocity  distribution 
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3.2  Free  jet  flow 

Figure  9  shows  the  experimental  set-up  of  a  free  jet  issuing  from  a  nozzle  and  the  Michelson  optics  used.  The 
nozzle  exit  diameter  is  d  =  2  mm.  The  jet  is  produced  by  blowing  carbon  dioxide  gas  out  of  a  pressure  chamber 
into  the  atmosphere  at  ambient  pressure.  The  tracer  particles  are  condensed  carbon  dioxide  particles. 


pressure  chamber 


Fig.  9  Schematic  of  the  nozzle  flow  and  optical  set  up 

The  beam  of  a  1  W  single-mode  Ar+-laser  is  transformed  into  a  light  sheet  with  a  cross-section  of  0.5  mm  x 
15  mm,  passing  axially  through  the  longitudinal  axis  of  the  jet.  The  angle  between  the  interferometer  and  the  flow 
direction  is  adjusted  to  90°.  A  glass  block  with  a  length  of  180  mm  and  a  refractive  index  of  n  =  1.52  is  used. 
Figure  10a  represents  a  Doppler  picture  of  the  jet  taken  with  the  CCD2  camera  at  a  pressure  in  the  pressure 
chamber  of  p  =  20  bar.  The  unshifted  reference  fringe  system  outside  of  the  jet  is  simultaneously  seen  in  Figure 
10b.  The  light  intensity  picture  in  Figure  11,  taken  by  the  camera  CCD1,  shows  the  Mach  disc  generated 
downstream  by  underexpansion.  Unfortunately,  due  to  the  shock  strength,  the  temperature  increases  and  the 
condensed  carbon  dioxide  particles  evaporate  in  a  small  region  behind  the  shock  wave.  Figure  12  shows  the 
Doppler  picture  (Figure  10a),  normalized  according  to  Equation  (2)  with  the  intensity  I0  image  ( Figure  1 1). 


Fig.  10a  Free  jet  Doppler  picture  from  CCD2  camera 


Fig.  1 1  Light  intensity  picture  taken  from  CCD1  camera 


Fig.  10b  Doppler  picture  of  the  jet  with 
reference  fringe  distribution  outside 


Fig.  12  Normalized  Doppler  picture  shown  in  pseudo  colors 

In  order  to  determine  the  velocities  (Figure  13),  the  technique  described  in  section  3.1  based  on  the  pixel  by 
pixel  gray  value  evaluation  is  carried  out  by  applying  Equations  (8)  -  (12).  Two  regions  had  been  evaluated:  a) 
shortly  after  the  nozzle  exit  and  b)  downstream  the  recompression  shock.  The  velocities  determined  by  applying 
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Equations  (12a,  b)  and  Equation  (1)  along  the  center  line  of  the  left  hand  side  region  a)  show  some  differences 
which  are  of  the  order  of  less  than  5%  compared  to  the  one-dimensional  theory,  see  Figure  14.  The  calculations 
give  a  higher  expansion  than  the  Doppler  picture  evaluation.  In  the  calculation  scheme  some  uncertainties  are 
included,  e.g.,  as  its  application  to  the  three  dimensional  jet  formation  as  well  as  the  not  exactly  known  pressure 
inside  of  the  pressure  chamber,  so  that  the  deviation  between  theory  and  the  velocity  distribution  taken  from  the 
Doppler  picture  is  acceptable.  Moreover,  the  Doppler  picture  shows  the  formation  of  the  radial  structure  of  the 
jet  in  the  domains  evaluated  and  visualized  in  Figure  13.  This  outcomes  support  the  good  applicability  of  the 
DPV  technique  for  planar  flow  velocity  determination  and  visualization. 


x  (mm) 

Fig.  14  velocity  distribution  along  the  jet  center  line 
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4.  High  speed  wedge  flow  visualization 


4.1  Description  of  shock  tunnel  facility 

The  DPV  technique  was  applied  to  measure  the  vertical  velocity 
component  over  a  wedge  in  supersonic  high  speed  flow.  That  flow 
is  generated  in  the  ISL/STB  shock  tunnel  facility  equipped  with  a 
conical  nozzle  used  in  direct  mode  and  producing  a  Mach  number 
4.0  nitrogen  flow  out  of  the  nozzle  exit.  Shortly  after  the  nozzle  exit, 
50  mm  downstream,  a  15°  wedge  is  placed  slightly  off  axis  (+7  ram) 
into  the  supersonic  nozzle  flow  producing  an  attached  oblique  shock 
wave  over  the  wedge  surface.  That  shock  wave  turns  the  incoming 
flow  parallel  to  the  wedge  surface.  A  principle  sketch  of  the  shock 
tunnel  facility  ISL/STB  in  which  these  experiments  have  been 
carried  out  is  given  in  Figure  15.  The  shock  tunnel  consists  of  three 
parts:  the  high  pressure  driver  section  with  driver  pressure  p5,  the 
driven  section  containing  the  test  gas  initially  under  pressure  p1;  and 
the  nozzle  attached  at  the  end  of  the  driven  tube.  As  driver  gas 
hydrogen  is  used  under  a  pressure  of  about  p5  «  130  bars.  The 
driven  tube  nitrogen  fill  pressure  was  p,  =  2,2  bars.  After  the  burst 
of  the  diaphragm  a  shock  wave  is  formed  by  which  the  test  gas  flow 
is  initiated  in  the  driven  section  and  expanded  via  the  shock  tunnel 
nozzle  in  order  to  form  the  wedge  on  flow. 

The  wedge  is  mounted  inside  a 
dump  tank  which  is  equipped  with 
windows  at  both  sides  and  at  the 
top.  A  light  sheet  is  produced  from 
top  by  a  Neodym/YAG  laser  pulse 
(523  nm),  see  vector  eL  in  Figure 
16,  which  illuminates  a  plane 
axially  oriented  in  the  middle  of  the 
wedge  flow.  The  observation  by  the 
Michelson  interferometer  was  done 
from  the  side,  eD,  and  the  flow 
vector  v  is  directed  vertically. 


Fig.  16  Vector  diagram 


Figure  17  gives  a  view  of  the  ISL  shock  tube  laboratory.  The 
shock  tube  used  for  the  experiments  reported  herein  is  that  on  the 
right  hand  side,  STB,  and  is  operated  as  a  shock  tunnel  in  non- 
reflected  mode.  On  the  left  hand,  shock  tube  STA  is  located.  The 
conical  nozzle  used  at  STB  opens  into  the  dump  tank,  see  Figure 
18,  with  the  wedge  mounted  in  front  of  the  nozzle. 


Fig.  1 7  Shock  tube  laboratory 


Fig.  18  Shock  tunnel  nozzle,  wedge 
arrangement  and  illumination 
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4.2  Wedge  flow  visualization  by  DPV 

In  Figure  19  the  initially  adjusted  fringe  pattern  present  before  flow  onset  can  be  seen  in  the  light  sheet  plane 
formed  over  the  wedge.  In  order  to  scatter  the  light  coming  from  the  top  into  the  Michelson  interferometer  in  case 
of  no  flow  (no  tracers),  a  diffuser  plate  was  placed  over  the  wedge  surface.  The  interference  fringe  adjustment  in 
tinning  mirror  Sj  is  clearly  visible.  The  fringes  appear  slightly  curved  due  to  optics  inaccuracies. 


Because  of  the  characteristics  of  the  conical  nozzle  used  for  the  experiments,  the  flow  around  the  wedge  is 
slightly  divergent  resulting  in  a  vertical  flow  component,  superimposed  to  the  axial  nozzle  flow,  which  is  in  the 
range  of  maximum  5%  of  the  axial  velocity.  Therefore,  the  fringe  pattern  in  front  of  the  oblique  wedge  shock 
wave  is  Doppler  shifted  as  well  as  behind  it  in  the  region  over  the  wedge  surface.  The  fringe  pattern  as  it  looks 
for  the  wedge  flow  is  given  in  Figure  20  (CCD2  camera).  The  wedge  shock  wave  is  clearly  visible  on  the 
Doppler  picture  by  means  of  a  fringe  shift  caused  by  the  velocity  jump  across  the  shock  wave  attached  to  the  tip 
of  the  wedge.  Over  the  wedge  the  flow  turns  vertically  by  flow  deflection.  That  vertical  velocity  causes  a  Doppler 
shift  of  the  frequency  of  the  light  scattered  by  tracer  particles,  here  titanium  dioxide  (Ti02)  which  results  in  a 
light  intensity  change  across  the  wedge  shock  wave  that  appears  as  approximately  one  half  fringe  displacement  in 
Figure  20.  The  processed  flow  region  is  marked  in  Figures  20  and  21  by  a  rectangle. 


Fig.  19:  Doppler  picture  taken  before  flow  onset 


Fig.  20  Doppler  picture  in  S’ 


As  already  seen  in  section  3.2,  in  gas  flows  the  tracer  particles  are  commonly  not  distributed  homogeneously 
resulting  in  a  distortion  of  the  Doppler  picture  light  intensity  distribution  I  in  the  image  plane  X'.  In  order  to  solve 
this  problem  independently  of  light  intensity  fluctuations  which  influence  the  amplitude  I0  in  Equation  (2),  a 
picture  of  the  light  intensity  I0  distribution,  see  Figure  21,  is  taken  simultaneously  with  the  Doppler  picture.  The 
optical  set-up  used  is  that  in  Figure  5.  In  dividing  the  two  pictures  pixel-wise,  the  Doppler  picture  light  intensity 
distribution  I  (Figure  20)  by  the  intensity  picture  I0  of  Figure  21,  a  normalized,  1 0  -independent  Doppler  picture 

results  according  to  Equation  (2)  with  0  <  I/I0  <  1.  This  intensity  normalizing  procedure  must  also  be  done  for  the 
Doppler  picture  taken  initially  before  flow  onset  (Figure  19). 


Fig.  21 :  Light  intensity  picture  in  X’  on  the  CCD1  camera 

At  each  pixel  of  the  CCD2  camera  (optical  set-up  see  Figure  5),  used  for  taking  the  Doppler  and  reference 
picture  the  light  intensity  is  stored  and  the  procedure  explained  in  section  3.1,  Equations  (8)  and  (9),  is  applied. 
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Then,  for  calculating  the  phase  difference  A(p,  on  the  one  hand  before  flow  onset,  Acpi,,  and  on  the  other  hand 
AtpD>  given  by  Doppler  shifted  scattered  light,  Equations  (10)  and  (11)  are  used.  By  subtracting  both  phase 
differences  from  each  other,  Equation  (12b),  the  frequency  shift  dvL  =  Vd  -  Vl  is  obtained  by  applying  Equation 
(12a).  For  calculating  the  amount  of  the  tracer  velocity  v  visualized  with  the  Doppler  picture  in  Figure  20,  the 
relation  deduced  from  Equation  (1), 

v  _  dv  l  c  1 

vL  2  cosy-sin 0/2  ’ 
is  applied  pixel  by  pixel. 

4.3  DPV  Velocity  picture 

The  image  processing  software  developed  at  ISL  applies  pixel-wise  the  procedure  described  in  section  3.1. 
Thereby  the  distribution  of  the  velocity  v  appearing  in  the  light  sheet  plane  was  determined.  The  result  is  given 
in  Figure  22  by  means  of  the  velocity  picture  which  shows  in  pseudo  colours  the  visualized  formation  of  the 
vertical  v-velocity  component  present  in  the  nozzle  flow  in  front  of  the  oblique  wedge  shock  wave  and  behind  it 
over  the  surface  of  the  wedge. 

Looking  at  the  velocity  evaluation  in  Figure  22,  deduced  from  the  Doppler  picture  in  Figure  20,  the  relative 
velocity  changes  Av  are  given.  Relative  to  zero  velocity  (v  =  0)  before  flow  onset,  the  vertical  nozzle  flow 
velocity  is  determined  to  Av  =  100  m/s.  The  measuring  error  is  in  the  range  of  about  10%,  i.e.,  ±  10  m/s.  That 
flow  region  appears  in  Figure  22  in  green  colour.  Across  the  wedge  oblique  shock  wave,  the  measured  velocity 
jump  is  of  the  order  of  Av  =  260  m/s  ±  26  m/s,  as  visualized  in  the  velocity  picture  of  Figure  22  by  a  red  coloured 
region  and  is  constant  between  shock  wave  and  wedge  surface. 


Fig.  22  Doppler  picture  visualized  velocity  changes  in  pseudo  colours 

4.4  Determination  of  nozzle  flow 

According  to  Equation  (10),  the  reference  fringe  pattern  in  Figure  19  gives  information  on  the  phase 
difference  distribution  Aqv  before  flow  onset.  Unfortunately,  probably  due  to  vibrations  present  during  the 
experimental  run,  the  reference  fringe  system  in  Figure  19  was  found  to  be  unstable.  The  reason  is  that  some 
optical  parts,  especially  the  mirrors,  are  displaced  resulting  in  a  shift  of  the  optical  path  length  difference  A<>,  see 
Equation  (3).  This  shift  is  much  smaller  than  one  fringe,  but  gives  a  uncertainty  in  the  Doppler  picture  evaluation 
by  applying  Equation  (12b)  concerning  the  phase  difference  Aq^.  It  was  found  that  the  Acp,  deviation  is  of  the 
order  of  one  quarter  of  a  fringe,  so  that  the  reference  picture  in  Figure  19  cannot  be  used  to  correctly  determine 
the  absolute  v-velocities  in  front  of  the  oblique  wedge  shock  wave  and  behind,  only  the  relative  ones.  To  obtain 
reliable  information  on  the  inexactness  of  the  DPV  determined  flow  velocity,  another  velocimetry  system  was 
applied,  the  so  called  Particle  Image  Velocimetry  (PIV).  Figure  23  shows  the  PIV-image  for  the  u-velocity 
distribution  in  horizontal  direction  and  Figure  24  that  for  the  vertical  v-speed.  A  detailed  description  of  the 
application  of  the  PIV  method  to  shock  tube  high  speed  flows  is  given  by  Havermann  et  al.  (2001). 
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In  Figures  23  and  24  the  15°  wedge  ramp  is  located  in  the  lower  part  of  the  diagram,  seen  in  white  colour. 
The  velocity  distribution  which  forms  in  the  light  sheet  plane  over  the  wedge  is  well  established  and  is  shown  in 
terms  of  pseudo  colours.  The  location  of  the  oblique  shock  wave  in  front  of  the  wedge  surface  can  be  determined 
as  a  strong  velocity  jump  in  both  pictures,  in  the  u-velocity  of  Figure  23  as  well  as  the  v-velocity  in  Figure  24. 
From  Figure  23  the  longitudinal  u-velocity  is  determined  in  front  and  behind  the  oblique  shock  wave.  Both 
values  are  deduced  as  mean  samples  coming  from  the  interrogation  regions  used  for  PIV  data  reduction.  Just  as 
done  for  the  axial  u-velocity,  also  for  the  vertical  v-velocity  the  data  are  deduced  for  the  onflow  and  for  the  flow 
over  the  wedge  surface.  The  following  data  are  obtained: 

front  (green/red  colours):  u  =  1257  ±  36  m/s,  behind  (blue  colours):  u  =  1 135  ±  33  m/s, 
front  (blue  colours):  v  =  35  ±  22  m/s,  behind  (yellow/red  colours):  v  =  305  ±  48  m/s. 


In  theory,  by  a  solution  of  the  conservation  equations  it  follows  that  the  shock  wave  influenced  wedge  flow  is 
constant.  Experimentally  this  prediction  is  confirmed  by  DPV  as  well  as  the  PIV  measurements,  see  Figures  22, 
23  and  24.  Moreover,  the  theoretical  prediction  gives  v  =  296  m/s  which  agrees  well  with  the  PIV  measurement 
(v  =  305  ±  48  m/s).  In  applying  ISL’s  shock  tube  flow  calculation  procedures,  the  flow  at  the  exit  of  the  nozzle 
can  be  determined  using  the  experimentally  obtained  data.  From  these  deductions  the  wedge  on  flow  is 
characterized  as: 

u  =  1257  m/s,  v  =  35  m/s,  M  =  4.0,  p  =  42  kPa,  T  =  240  K. 


Fig.  24  Vertical  velocity  distribution 
measured  with  PIV 


x/ mm 
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4.5  DPV  reference  picture  calibration 

As  discussed  in  section  4.4,  the  PIV  velocity  visualization  in  Figure  23  shows  that  the  wedge  onflow 
produced  by  the  conical  shock  tunnel  nozzle  flow  is  superimposed  by  a  vertical  velocity  component  of  v  =  35  ± 
22  m/s.  Looking  at  the  velocity  evaluation  in  Figure  22  deduced  from  the  Doppler  picture  in  Figure  20,  there  the 
vertical  nozzle  flow  velocity  was  determined  to  v  =  100  ±  10  m/s.  This  outcome  does  not  agree  with  the  v- 
velocity,  v  =  35  m/s,  measured  with  the  PIV  technology  which  is  successfully  tested  by  Havermann  et  al  (2001). 
The  reason  for  the  differences  in  the  Doppler  picture  evaluation  is  already  discussed  in  section  4.4  as  a  result  of 
optical  length  A<>  variations  between  taking  the  reference  Doppler  picture  (Figure  19)  and  Doppler  picture 
(Figure  20),  due  to  small  vibrations  of  the  optical  parts  inside  the  Michelson  interferometer.  In  consequence,  a 
systematic  error  develops  which  must  be  taken  into  account  as  a  corrective  speed  factor.  For  the  experiment 
described  in  section  4.2,  that  factor  results  to  65  m/s.  In  subtracting  this  corrective  factor  from  the  velocity  data 
given  in  Figure  22,  the  absolute  velocity  distribution  can  be  obtained  as  shown  in  Figure  25. 


Fig.  25  Calibrated  DPV  velocity  picture  in  terms  of  pseudo  colours 

If  the  reference  interference  picture  is  taken  simultaneously  with  the  Doppler  picture,  as  done  in  the 
experiments  described  in  sections  3.1  and  3.2,  the  absolute  velocity  distribution  is  available  with  a  measuring 
error  of  about  10%.  Due  to  mechanical  instabilities,  in  case  the  reference  picture  is  taken  before  taking  the 
Doppler  picture  with  the  optical  set-up  shown  in  Figure  5,  only  velocity  variations  are  given  correctly,  as  the 
velocity  jump  across  the  wedge  shock  wave  in  the  experiment  carried  out  in  the  ISL  shock  tunnel  STB.  This 
jump  is  measured  to  Av  =  260  ±  26  m/s.  The  PIV  data  in  Figure  24  give:  Av  =  270  ±  35  m/s.  This  outcome 
suggests  that  the  DPV  technique  at  moment  is  only  able  to  measure  correctly  velocity  changes  (e.g.  the  jump 
across  the  oblique  wedge  shock),  because  the  reference  and  Doppler  pictures  are  not  simultaneously  available. 


5.  Conclusions 

The  Doppler  Picture  Velocimetry  (DPV)  technique  began  in  1981  with  the  visualization  of  fringe  shifts  showing 
changes  of  flow  velocity.  At  that  time  the  fringe  shift  determination  only  gave  information  on  the  velocity 
changes  in  the  flow  area  considered.  Nowadays,  in  the  enhanced  DPV  technique,  the  Doppler  picture  is  taken 
with  a  CCD-camera  and  the  output  is  treated  pixel-wise  using  an  appropriate  software  which  is  under 
development  at  ISL.  The  DPV  concept,  using  the  gray  values,  is  a  very  powerful  tool  allowing  to  obtain 
information  at  any  pixel  of  the  plane  crossing  the  flow. 

To  overcome  the  uncertainties  concerning  the  reference  picture  which  is  taken  before  flow  onset  and 
discussed  in  section  4.2,  a  new  Michelson  system  is  established  for  taking  the  Doppler  picture  (phase  difference 
A<Pd)  and  the  reference  picture  (phase  difference  A(pL)  simultaneously.  For  that  purpose,  the  optical  set-up  given 
in  Figure  26,  is  under  development  and  will  be  tested  in  future  applications  of  flow  velocity  measurements  in  the 
ISL  shock  tunnel  facilities.  The  new  DPV  tool  uses  three  CCD  cameras,  CCD1  for  taking  the  interference  fringes 
resulting  from  the  Doppler  shifted  frequency  of  the  scattered  light,  CCD2  for  taking  the  reference  picture  which 
uses  unshifted  laser  light  and  CCD3  for  taking  the  intensity  distribution  of  the  scattered  light.  The  light  passing 
the  Michelson  interferometer  is  polarized  in  order  to  separate  the  Doppler  shifted  light  and  the  reference  light 
which  comes  from  a  dispersion  plate  that  is  directly  illuminated  by  the  laser  light  source. 
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Fig.  26  Optical  set-up  for  simultaneously  taking  Doppler  (CCD1)  and  reference  picture  (CCD2) 
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